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Abstract 

Starting from the parametric representation of a Feynman diagram, 
we obtain it's well defined value in dimensional regularisation by chang- 
ing the integrals over parameters into contour integrals. That way we 
eventually arrive at a representation consisting of well-defined com- 
pact integrals. The result is a simple transformation of the integrand 
which gives the analytic continuation of a wide class of Feynmaninte- 
grals. The algorithm will especially be fit for numerical calculation of 
general massive multi-loop integrals. An important advantage of this 
method is that it allows us to calculate both infinite and finite parts 
independently. 
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1 Introduction 



In recent years there has been an increasing interest for the evaluation of 
massive multiloop Feynman diagrams. High precision experiments force the 
theoretical predictions to reach an equal amount of accuracy. Over the years 
several methods have been proposed to deal with the problem. 

The first succesful attempts to calculate Feynman-diagrams systemati- 
cally were integration by parts |l| , IR-rearrangement method and gegen- 
bauer technique |Q. They are not fit though for calculating finite parts of 
general massive diagrams, although they have led to some quite impressive 
calculations of beta- functions in different theories Q . An excellent review 
of these methods is given in Q. 

A first general massive approach |7[ ^ ^ is analytic by nature and 
is based on the following basic principle: by putting a certain number of 
masses in the diagram equal to zero a gamma-function ansatz in obtained, 
upon which the masses can be added again by means of the Mellin-Barnes 
representation. This gives rise to hypergeometric series which have proven to 
be rather succesful, e.g. for the asymptotic expansions in the two loop case 



[10|. These series have certain drawbacks however: they converge only in 
certain kinematic regions and although the ansatz needed for the application 
of this method can easily be found in the case of two loops, more loops will 
give considerable problems. Applications of these methods are therefore 
mainly confined to asymptotic 2-loops cases. 

Other algorithms avoid the use of Euler-gamma functions and there- 
fore turn out to be numerical. An excellent example of such a method, 
which even avoids parametric representation and Wick-rotation, relies on 
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the seperation of an orthogonal space of momenta and integrating out this 
space first. Although this method has been succesfully used for some spe- 



cific cases [11, 12, 13, 14 1, it has not yet been expanded up to three loops 



or more. Other numerical aproaches [15, ITH[ relies on the parametric repre 



sentation of the diagram. These approaches also differ mainly form ours in 
this respect that these methods are based on a subtraction of divergences 
under the parametric integral, while ours is fundamently based on dimen- 
sional regularisation. Analytic continuation in d dimensions is central in our 
approach: thus we obtain not only the finite parts but also the poles of the 
Feynman-diagram. It is therefore that it has no trouble in dealing with IR 
and UV divergences at the same time. Our method is without reservation 
applicable in a vast number of cases and easily implemented on a computer. 

The rest of the article is organised as follows. Section two contains the 
basic formulas of the contour method and some adaptations desirable for 
smooth numerical calculations. Section three consists of some examples and 
section four is a summary and conclusion. 

2 Contour method 

A scalar diagram D with L loops, / internal lines each labeled by number / 
and a mass mi, a total external momentum Pj per vertex and a space-time 
dimension d has a parametric representation [^^ 

lD{P)=i rida, [(4,p.^^(,)]./2 (1) 







where 



QD{P,a) = . X! n (3) 

with T the set of all the trees of D and the set of all the two-trees of 
D and the square of the momentum which passes through the cut. 
is a homogeneous polynomial in a of degree L and of degree 1. If the 
diagram gives rise to tensor integrals with irreducible numerators we can 
change these to scalar integrals using for example the general expression in 



mm- 

The general philosophy of our method will be as follows: we will isolate 
the different poles in the integrand and then avoid them by changing the 
integral in a contour integral flung around the pole. This way we will obtain 
a well-defined analytic continuation of the Feynman integral in dimensional 
regularisation. 

Our first goal will be to isolate the poles in (|l|). Since the polynomial 
Rd of e.g. the setting-sun diagram has the form 

this is not allways obvious: no simple poles are visible. There is however a 



well-known substitution, used in the convergence theorem |17|, which does 



the job. It involves a separation of the integration domain into sectors 
< a<^(i) < a^(2) < • • • < a„(i) 

with a a permutation of (1, 2, . . . , /). Per sector we preform the following 
change of variables 

"tT{i) = • • • I32P1 

a<7(2) = PiPi-i ...(32 



a^ii) = Pi (4) 

where 

< /?/ < oo 
0<A<1 / = !,.. .,1-1 

Now we can interpret each sector as a family of nested subsets of Unes of 
the diagram D: these nested subsets Si C Sj are determined by the scaHng 
behaviour of (3i and f3j, i > j. 



We can prove the fohowing form for Ru [17| 

Rom = (3^'f3t . ■ ■ pf'ii + om (5) 

where Li is the number of independent loops in the subset of lines spec- 
ified by the lines which disappear if the corresponding Pi is put to zero 
(remember: each line corresponds to a certain a). Algebraically we have 
independent poles in the denominator and now we will be able to preform 
our " contouration" . 

The parametric representation (||) becomes 

exp -E. UU + QUP^ /^)] 

thus a sum over permutations cr, which means a maximum of /! terms. 
Thanks to the symmetry in the diagram a certain collecting of terms will 
usually be possible. From now on we will omit the sommation over a and 
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concentrate on one sector only. For notational convenience we take a to be 



the identical permutation. Due to homogeneity of Rd and Qd [17| we can 
immediately preform integration over /?/ and get 



(47r) 



(6) 



[i?i,(/3i,...,/3,_i)]'^/2 
We will omit these independent factors in front of the integration in what 
follows. 

Only Ry) gives rise to possible poles if the diagram is globally UV- 
divergent (i.e. L(i/2 — / > 0). If the diagram is UV-convergent and some 
masses are zero then poles can come from the numerator in (which is then 
in the denominator!): they correspond to IR-divergences. If L(i/2 — I > 
and masses are zero, these IR-factors may compensate some UV-poles: this 
is a well known fact in dimensional regularisation. Now we can write Id in 
the following form 

Id{P)= f'dpi--- TdA-i f'dPi+i--- f'dPi-i 
Jo Jo Jo Jo 

d/3i • • • / dA-1 / dA+1 • • • / d/?/_i 



JO JO JO 

■hi) (7) 

Pi " ' Pl-l pl-1 ' ' ' PI-1 
where the function /(/3i, . . . , (3i) is written as as /(/3) for convenience. We 
define a unique pi and / by demanding that /(/?» = 0) 7^ and that it is 
analytic at /3j = 0. In dimensional regularisation, i.e. d = 4 — 2e, 

Pi = ni- QiE, 
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c 



1 



Figure 1: The contour for obtaining analytic continuation of /3-integrals 

rii integer and qi rational. From now on we will omit the cases nj < because 
in these cases the integration can be preformed without regularisation. 

Let's concentrate on the integration of one specific (3i. In the function 
/(/?) all other (3^, i ^ I will be considered to be parameters for the time 
being. We will use the notation f{Pi)- The function /(/?;) is analytic at the 
origin by definition. In dimensional regularization (i.e. qi ^ 0) the integrand 
as a whole is analytic in a region around the real axis for Re/?; > minus 
the cut line along this positive axis. This is true for every / = 1, . . . , I. 

We will obtain a meaningful regularized value for the integral if 
we change the integral into a contour-integral, which conincides with the 
original integral in non-divergent cases. In concreto we define the Ui (un- 
renormalized) operator (e.g. |l^]) to be 

where the contour C, as is showed in figure 1 , indeed avoids the pole A = 
by a tiny contour Cq which is a circle with radius A, limA — > 0+ and the 
integrand is analytic in the region of the contour as shown before. It is clear 
that this definition only holds for non- integer pi: our method is intimately 
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connected to dimensional regularization. 

In order to get rid of contours again, we preform integration over Co 
explicitly. We are abel to do so if we write /(A) in a Taylor expansion (this 
is allowed due to analyticity of / at the origin). Since limA 0+ we keep 
the ni first terms of the Taylor expansion. Other terms vanish in this limit, 
they are namely 0(A^'^) and higher. Finally we obtain 



This form still contains A explicitly. If we use the identity for every pi ^ 



Noticing that the integrand under the first integral as a whole is of the order 
f3i'^ , the integral will be convergent in the limit A ^ 0. We will use the 
following notation 





and regroup the terms, we find for the case rii > 1 





and for the other terms 




j! \j -Pl + 1 



) 



zi=0 



They are independent of (3i. 
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We can write this result as a simple transformation of the integral 



'dA/f^+E/f? (9) 
j 

The next difficulty we have to face, is whether or not this operation can be 
repeated for other /?'s? Potential problems arise when + Qd{P,(3) 

comes in the denominator: we have not been able to prove an expression 
like (^) for it, but this appears to be always the case. We have done numer- 
ous tests on diagrams up to three loops under different circumstances (sev- 
eral external momenta, masses or massless,...) and have found no counter- 
examples. Yet a rigorous mathematical proof that this applies for every 
diagram is still lacking. Here we will assume that the diagrams we are 
dealing cause no problems. 

In that case if we have carried out the C/^j -operation, we can re-establish 
the other /3j as true variables. Then we choose another I2 to focus on, write 

1(1, u\ — I dpu n. 



and apply the operator Ui^ to I(i^^i^y We can maintain our notation by 
adding a number to the sub- and superscript of /, e.g. after U2 o Ui we will 
get 

f _ f (*:*) I f («.*) 1 f (*j) I f (* j) 

/ (h,l2) - J {h,l2) ^ 2^J ihM) ^ {h,h) ^ {h,l2) 

We repeat this procedure until no /?-poles are left, thus regulazing the 
feynman- integral Id{P) completely. 
Several remarks are in order. 
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• We can regulaze every variable /?» before preforming integration, i.e. 

commutes with / d/3j for « 7^ j. We only have to take care to 
keep the different parts of /j-^^ (all other indices are arbitrary: j or 
*) together, because each term on his own has a pole in Pi, only the 
sum converges. Before preforming integration we can also expand in 
e. This allows us to calculate the different coefficients of the laurent- 
expansion in e seperately. Thus the T operation of [l^ (isolating the 
divergent part) is easily implemented. Higher pole parts can only come 
from the f^j^ parts. 

• The result of the different operations C/;^ °Ui2°- ■ is independent of 
their order. If there are only poles in f5i and /32, so Ui^ °Ui^ completely 
regularizes the integral then: 

This follows from the fact that the [/-operators are in fact analytic 
continuations (see our primary definition (^)): ifHepi < 1 the contour- 
integrals coincide with the ordinary integrals which are convergent, so 
Ui-^ o = Ui^ o Ui-^ in the area Repi < 1. By the principle of analytic 
continuation they must coincide for all complex d (which fixes all the 
Pi's), for which Ui-^ and Ui^ exist. Of course this implies that for every 
variable (3i this analytic continuation has been carried through. 

It is clear that an expression like (^) will be especially useful in numer- 
ical calculations and therefore it might be useful to examen it's numerical 
behaviour. Only /^^-^ will cause trouble around zero. Although the func- 
tion as a whole is finite, it is really a subtraction of one or more infinite 
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values. Moreover it is of the order /3;^'^ which after expanding in e gives 

1 + qieln{(3i) -\ and is thus another thread to numerical stability. There- 

fore we will approximate f^^^ around zero, by the next term in the Taylor 
expansion, i.e. 



zi=o nil 



'(0 ~ d-^^zi 

Introducing this approximation in the interval [0, (^], 5 -C 1 being a certain 
numerical value, we obtain 



"'/3f' d^izi 



^I'^-W (10) 



The intrinsic error of this approximation depends on the highest order pole. 
If it is l/e^, then terms proportional to ln{5)^^^ will appear. If 5 is very 
small then this correction will be large, resulting in a smaller accuracy. Thus 
there is a best value for S typical for every diagram. Typical values are e.g. 
6 = 10~^ . In order to obtain larger values for 6 without losing accuracy one 
could include more terms of the taylor expansion. 



3 Examples 

In this section we will discuss one example at length: the famous "setting- 
sun" diagram in scalar A^^. Then some other results for diagrams with a 
larger number of propagators are included. In order to be able to compare 
with analytic results, we begin by choosing some trivial vacuumbubbles: no 
external momenta and equal masses. Of course this is only to make compar- 
ison possible: it is the analytic methods that are bound by these limitations, 
not the contour method as will be shown in some explicit numerical exam- 
ples. 
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We put the momentum p equal to zero and choose the masses equal. This 
makes the diagram is very symmetrical: topologically every leg is equivalent 
and this will simplifie the explicit calculation. 

Starting from the well-known parametric representation of the setting- 
sun diagram, we get only one /3-sector (all masses are equal). Applying 
formula (^)), we get 

,2\d-3 rl rl ^ ^ R.\d-3 



'ss- 



-T{3-d)3\ d(3i ^ , ^, 11 



The diagram is globally divergent in d > 3 whence the factor T{3 — d). At 
this stage we see a simple pole in /?2 emerging. This corresponds to the 
two-legged subdivergence present in the diagram (renormalization of the 
coupling constant). With the notations of section 2, we have 

and P2 = n2 — q2£ = 1 — e. Applying formula (^) we get 
{l+(32 + (3il32Y-'' 1 \ 



TT f (1+/?2 + /?1, 



L 



/'(l + /32 + /3i/32)^-=^ 

<iP2 



+- 



(1 + A + PiP2r/' (1 + Pi)"/^ J p^/^-' 
1 



(2-d/2)(l + /3i)'^/2 
These integrals are finite, as expected. If we restrict ourselves to the diver- 
gent part (without expanding the r(3 — d) though, thus fj^^) ^^'^ the fj^2)^s 
up to e") we get 

/« = [\P2^ ^^P^ + ^^h 



(2) io "\il + Pl+(3lP2? {l+HlYJPl 

M) _ 1 , ln(l + /3i) 

J{2) ~ , a \9 "T 



e(l + /3i)2 {l + Pif 
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We can preform this integration analytically 




If we expand further in e and preform the numerical calculations we obtain 
the finite part as well: 

C^2\d-3 /o \ 

/gg = \ > r(3 - d) - + 3 - 8.966523919e (12) 
(47r)'* \£ J 

Comparing this result by numerical evaluation of the analytic results in ||^] 
we see that only the last digit is different. 

Another test we did was to repeat some of the cases from i.e. a sub- 
tracted setting-sun diagram with arbitrary masses and external momentum. 
A selection 



p2 


mi 


1712 


1713 


Contour 


Series 


9 


3 


3 


10 


-7.31299 


-7.31298 


49 


1 


1 


10 


-0.316742 


-0.31675 


-25 


2 


2 


10 


-1.942844 


-1.94285 


-250 


4 


4 


4 


-13.57190 


-13.5719 


100 


3 


3 


3 


-1.282852-i20. 89960 


-1.28285-i 20.8996 


49 


20 


20 


10 


-1014.695748 





The last result cannot be calculated with the method of Q because mi -|- 
m2 < 711-3 should apply under the threshold. 
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Another test for our method was the so-called "basketball-diagram", a 
vacuiimdiagram in with 3 loops and 4 legs. We also assumed equal 
masses here in order to be able to compare to known analytic results. We 
have calculated the diagram in four cases: with one, two, three and four 
massive legs. 

The numerical procedure as given by formula (|l^) in section two, was im- 
plemented in MATHEMATICA the diagram, expressed in it's original 
form is automatically transformed into the parametric o-form and /3-form, 
then cast into the expression ( |lO|) and numerically evaluated. Because we 
are dealing with a very limited class of functions, it is possible to program a 
package that preforms one integration analytically. This general procedure 
was used to evaluate all the following results. 

These were the results we obtained up to order e^. We compare these 
values to analytic results in [20|. All results are given up to a factor 
(4^^3d/2 r(4 — 3d/2). Numerical evaluations of the analytic results are 
only shown in relevant cases. 



Type 


Contour 


Analytic 










el 


eO 




1 mass 




-0.5 


-1.5 


-5.967401126 


-1.5 


-5.967401100 


2 masses 


2 


-2 


-19.86960443 


22.33542372 


-19.86960440 


22.335425302 


3 masses 


6 


- 4.5 


-55.10881324 


44.04808956 


-55.10881320 


44.04808897 


4 masses 


12 


-8 


-107.2176265 


26.87588164 


-107.2176264 


26.87588031 



These results were obtained with 6 = 10 ^ and relative errors are of or- 
der 10-s. 

Other diagrams we evalutated are the 3-loop, 5 propagator diagram with 
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2 masses 

h = I dk ^ 



and the mercedes diagram with 3 masses in a very asymetric configuration 
h = fdk ^ 



'kl{kl + m?)kl{{ki - k^f + m'^){ki - k2Y{{k2 - k^f + m^) 

where dk stands for integration over ah internal momenta ki. As a result 
we found: 



(m2)5-3'^/2r(5 - 3d/2) \ 1 3 



(47r)3<i/2 

and 



2 + - + 2.065197729 + 19.42850404e 



(m2)6-3rf/2r(6 - 3d/2) , 

^ ' ^ '—L [7.21234140 - 9.081028480e] 



(4^)3d/2 

In the case of Is we encounter relative errors of order 10~^ in comparison with 



analytic results. According to ||2^ Iq has not yet been evaluated analytically. 
Our numerical results differ only by order 10^^ from theirs. We can conclude 
that the numerical behaviour of the method is fine. 



4 Summary and conclusion 

The main merit of our contour-method as introduced here, lies evidently in 
its general character. Indeed it is applicable for a wide class of Feynman- 
diagrams scalar and tensorial (although the borders of this class or not 
yet quite clear). The computer-algorithm based on expression ([lO| ) is easy 
to implement and virtually independent of the actual diagram fed to the 
system. 
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In the section examples we confined ourselves mainly to diagrams for 
which analytic expressions are known, in order to investigate numerical be- 
haviour. As shown in the case of the setting-sun diagram our contour- 
method is without further trouble applicable in the cases of external mo- 
menta and abitary masses and that is where her strength lies. 

As remarked our contour-method is a numerical method in the first place. 
The necessary calculations involve especially a lot of numerical integrations, 
which are known not to be an easy problem. Even at this point our method 
has some advantages: the integration intervals are compact at all time (^) 
and -thanks to our numerical adjustment (p!o|)- all integrands are finite. Such 
calculations can best be done by some adaptive Monte-Carlo method. The 
method will also give rise to a large number of integrals, but this drawback 



can partially be met using the recursion algoritms in |23, |2^, forcing 
us only to calculate a small number of so-called master-integrals. This is 
a common practice: such techniques have to be used in any method in- 
volving higher-order Feynman-diagrams in order to reduce the number of 
calculations. 

To sum up, we have a method which is easy the implement and applica- 
ble in a vast number of cases, even those where tradional methods fail, such 
as finite parts of complicated diagrams. Furthermore we are able to cal- 
culate the different coefficients of the laurent-expansion in e independently. 
Therefore it is a useful tool for high-precision calculation of general massive 
Feynman-diagrams. 
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